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Pure Yang-Mills SU(N) theory Is studied in four dimensional space and Landau gauge by a double 
perturbative expansion based on a massive free-particle propagator. By dimensional regularization, 
all diverging mass terms cancel exactly in the double expansion, without the need to include mass 
counterterms that would spoil the symmetry of the original Lagrangian. The emerging perturbation 
theory is safe in the infrared and shares the same behaviour of the standard perturbation theory in 
the UV. At one-loop, gluon and ghost propagators are found in excellent agreement with the data of 
lattice simulations and an infrared-safe running coupling is derived. A natural scale m = 0.6 — 0.7 
GeV is extracted from the data for N = 3. 
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Our knowledge of QCD is still limited by the lack of a 
powerful analytical method for the study of the infrared 
range. It is widely believed that perturbation theory 
breaks down at the low-energy scale Aqcd ~ 200 MeV 
and the study of very important phenomena, including 
hadronization and quark confinement, must rely on phe¬ 
nomenological models or numerical lattice simulations. 
Even pure Yang-Mills SU(N) theory is still not fully un¬ 
derstood in its infrared limit. Important progresses have 
been achieved in the last years by developing new an¬ 
alytical to olstl| 431| and by simulating larger and larger 
lattices |32H34l |. 

A key role is played by the dynamical massj^ that 
the gluon seems to acquire in the infrared according to al¬ 
most all non-perturbative studies. Moreover, it has been 
shown that the inclusion of a mass by hand leads to a 
phenomenolo gica l model that can be studied by pertur¬ 
bation theory [Hiii. 


In a recent paper [^, we pointed out that a mas¬ 
sive perturbation theory, based on a massive free-particle 
propagator can be developed by an unconventional set¬ 
ting of the perturbative method, without changing the 
Lagrangian and without adding free parameters that 
were not in the original Lagrangian, yielding a first- 
principle anlytical treatment based on perturbation the¬ 
ory. 


In this paper the idea is developed by a double ex¬ 
pansion in dimensional regularization. Pure Yang-Mills 
SU(N) theory is studied in four dimensional space and 
Landau gauge by an expansion in powers of the total 
interaction and of the coupling. By dimensional regu¬ 
larization, all diverging mass terms cancel exactly in the 
double expansion, without the need to include mass coun¬ 
terterms that would spoil the symmetry of the original 
Lagrangian. The emerging perturbation theory is safe in 
the infrared and shares the same behaviour of the stan¬ 
dard perturbation theory in the UV. While the present 
letter summarizes the main results, more details of the 
calculation will be published elsewhere [s^. 


Let us consider pure Yang-Mills SU{N) gauge theory 
without external fermions. The Lagrangian can be writ¬ 
ten as T = Lym + 'Clfix where Tvm is the standard 
Yang-Mills term and L fix is the gauge fixing term. Usu¬ 
ally the total Lagragian is split into two parts, a free- 
particle Lagrangian Lq that does not depend on the cou¬ 
pling strength g, and an interaction Lint that contains 
0(g) and O(g^) terms. In Landau gauge the gluon free- 
particle propagator is transverse 

A^‘^(p)=t'^^(p)^; = ( 1 ) 

and has a pole at p^ = 0. We can shift the pole and 
modify the free-particle propagator as 


Ai;-'(p) ^ , (2) 

—p^ -h 

without changing the content of the theory provided that 
the counterterm 


SL = (3) 

is added to the interaction. Actually, we are just adding 
and subtracting the same quantity in Tq s-nd Lint with¬ 
out changing the total Lagrangian. Thus we can develop 
a perturbative expansion in powers of the total inter¬ 
action and use the standard formalism of perturbation 
theory with a total interaction Lint + i^L that is a mix¬ 
ture of terms that depend on the coupling strength g and 
the counterterm that does not vanish in the limit g —>■ 0. 
In the expansion, the free-particle propagator is the mas¬ 
sive propagator ®, and the vertices are of order O(g^) 
(the counterterm), order 0(g) (the three-particle gluon 
and ghost-gluon vertex) and order 0(g'^) (the four-gluon 
vertex). That the content of the theory has not changed 
can be easily seen by summing up all graphs with n in¬ 
sertions of the counterterm in a gluon line. As shown in 
Fig. I, we can write a dressed propagator as the inhnite 
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Figure 1: Graphical illustration of Eq.(l4]). The cross is the 
counterterm of Eq. m that gives a factor m^. 


sum of a set of reducible graphs 
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The same result can be described as the effect of a proper 
polarization term If = rn? that arises from the couiitert- 
erm and cancels the mass. Then formally, the two ex¬ 
pansions are equivalent if we sum up all graphs. On the 
other hand, by the same token, at any finite order, the 
perturbation theory that we develop is not equivalent to 
the standard perturbation theory, but the two expansions 
differ by an infinite class of graphs that introduce non- 
perturbative effects. In other words, the two expansions 
differ by some non-perturbative content at any finite or¬ 
der. In fact, it is well known that the gauge invariance 
of the theory does not allow any shift of the pole in the 
propagator at any finite order, so that the massive zeroth 
order propagator Am cannot be obtained by the standard 
perturbation theory at any finite order. 

Since we already know, by non-perturbative calcula¬ 
tions, that the gluon propagator is massive in the in¬ 
frared, then we expect that the present expansion, with a 
massive zeroth-order propagator, would be more reliable 
than the standard expansion, at least in the infrared. 

While massive models have been studied before and 
found in fair agreement with the lattice datajs^-H^, the 
present approach is very different because the Lagrangian 
is not modified and does not break BRST symmetry, so 
that no free parameters are added to the exact Yang- 
Mills theory, yielding a description that is based on first 
principles and can be improved order by order. For in¬ 
stance, we can easily show that no mass is predicted for 
the photon by the same method. 

Having left the Lagrangian unchanged, we expect that 
the massive expansion should share the same behaviour 
of the standard expansion in the UV where any finite 
mass becomes negligible. In fact, if ^ the geomet¬ 
ric expansion in Eq. (|3]) is convergent and the two pertur¬ 
bation theories must give the same results. On the other 
hand, when p^ —>■ rn? each single term of the geometric 
expansion Eq.Q diverges and the formal sum of infi¬ 
nite poles introduces some non-perturbative content (^. 
Thus we can predict that the scale m should be close 
to the Landau pole A where the standard perturbation 
theory breaks down. 

An other interesting aspect of the present massive ex¬ 
pansion is that no other mass counterterms are required. 


Thus there is no need to include terms that, because of 
gauge invariance, were not in the original Lagrangian. 
All diverging mass terms are cancelled exactly by the 
counterterm 5L. We expected that a cancellation like 
that would arise by summing all graphs, just because the 
Lagrangian has not been changed and no diverging mass 
terms are present in the standard expansion. However, if 
we inspect the graphs in Eig.2, we can easily see that any 
insertion of 5L in a loop reduces the degree of divergence 
of the graph, and all mass terms become finite by a fi¬ 
nite number of insertions. Thus, if the divergences must 
cancel, they will cancel at a finite order of the expansion 
provided that we retain more counterterm insertions than 
loops. 

As pointed out before, the order of a graph is the num¬ 
ber of vertices that are included, while the number of 
loops is equal to the powers of in the graph. If the 
effective coupling is small, as it turns out to be accord¬ 
ing to non-perturbative calculations, we could consider a 
double expansion in powers of the total interaction and 
in powers of the coupling: we can expand up to the or¬ 
der n, retaining graphs with n vertices at most, and then 
neglect all graphs with more than ^ loops. If n is large 
enough, then all divergences in the mass terms are can¬ 
celled by the counterterms in the loops. Eor instance, at 
one loop we only need n = 3. 

In this paper we report the results for a third-order 
one-loop double expansion in dimensional regularization. 
The gluon polarization and the ghost self-energy are eval¬ 
uated by the sum of all graphs with no more than three 
vertices and no more than one loop, as shown in Eig.2. 
The integrals are evaluated analytically by dimensional 
regularization and expanded in powers of e = d — 4. 

In dimensional regularization, the cancellation of all 
diverging mass-terms can be easily proven by a simple 
argument. The insertion of just one counterterm in a 



(1a) (1b) (1c) (Id) 



Figure 2: Two-point graphs with no more than three vertices 
and no more than one loop. The ghost self energy and the 
gluon polarization contributing to the functions F and G are 
obtained by the sum of all the graphs in the figure. 
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loop can be seen as the replacement 


and adding the crossed loop with its symmetry factor 


1 1 2 1 2 d 

—+ rrfi -(- rn? ™ —p^ + m? ^ dvn? ™ 

(5) 

in the internal gluon line. If the graph has no other 
counterterm insertions, then its dependence on comes 
only from the massive propagators and a derivative of 
the whole nth-order £-loop graph gives the sum of all 
(n -I- l)-order .^-loop graphs that can be written by a sin¬ 
gle insertion in any possible ways. According to Eq.®, 
each diverging mass term rn^/e that comes from a loop 
would give a crossed-loop term —rn^/e. The argument 
also suggests a simple way to evaluate the crossed-loop 
graphs by Eq.®. 

The cancellation of all diverging mass terms without 
inclusion of any other conterterm is very important be¬ 
cause there is no need to include free parameters that 
were not in the Lagrangian, while all other divergences 
can be dealt with by standard wave function renormal¬ 
ization. 

It is instructive to inspect the constant graphs that 
contribute to the proper gluon polarization 11 in Fig 2. 
At the lowest order (n = 1, £ = 0) the counterterm 5L 
just adds the constant term Ilia = that cancels the 
shift of the pole in the propagator. The tadpole Ilif, is 


Ilib = -;Ct m 


2 , 1 
- + log — + a 

e D 


( 6 ) 


where the effective coupling a is given by 


SAT 

47r 


47r 


(7) 


The crossed tadpole IIic can be evaluated directly or by 
a derivative according to Eq.® 


Hie = -TO 


2 dlixb 3 2 

5to2 “ 4 “"" 



The diverging terms already cancel in the sum IIi;, + 
IIic. In fact, the double-crossed tadpole Hid is finite and 
including its symmetry factor it reads 

nid = -^aTO^ (9) 

O 


so that the sum is 


Ilib -f IIic -f Ilid — —a rn?. (10) 

o 

While the ghost loop vanishes in the limit p —>■ 0, a mass 
term can arise from the gluon loop that in the same limit 
is 


'2 

Il 2 b( 0 ) = -am^ ( - -I- log At + const. 

e TO^ 


( 11 ) 


n2b(0) + n2c(0) =(^1- n2b = -am^. (12) 

Thus the one-loop gluon propagator reads 

A(p)"^ = -p^ + ^arri^ - [n(p) - n(0)]. (13) 


We observe that a finite mass term has survived, and it is 
of order a. Since it only arises from the gluon loops, no 
mass would survive in QED for the photon by the same 
method. 

The calculation of the total one-loop polarization and 
of the ghost self energy is straightforward but tedious. 
For a massive theory, the one-loop sunrise graphs n 2 a (p), 
Il 2 b(p) and the one-loop ghost self energy have been eval¬ 
uated by several authors. The crossed loops follow by 
a mass derivative according to Eq.([S]). In the minimal 
subtraction scheme, all the divergences are cancelled by 
the same wave function renormalization constants of the 
standard expansion. Namely we obtain 


Za-1 = 


13a 


- 1 = 


a 


(14) 


The resulting propagators are finite, and the ghost and 
gluon dressing functions, x J respectively, can be 
written in units of the scale m?. Without taking any 
special subtraction point, the dressing functions can be 
recast as 


J{s) 


_’7(si)_ 

1 -I- aJ(si) [F{s) - E(si)] 


(15) 


X{s) 


_X(52)_ 

1 -f ax(s 2 ) [G(s) - G(s 2 )] 


(16) 


where s = p’^jm'^. The integration points si, S 2 are 
arbitrary as also are the normalization constants J(si), 
x(s 2 )- The functions F{x), G{x) do not depend on any 
scale or parameter and are given by the following explicit 
expressions 


5 1 

Fix) = — -f — \La + Lb + Lc + Ra + Rb + Rc\ 
G(x) = ^[Lg + Rg] (17) 


where the logarithmic functions Lx are 


^ , 3x3 - 34x2 - 28x - 24 

La[x) = - X 


4 -I- X / \/4 -f X — \/x 


■log 


\/47 


\/x 


Lb{x) = 


2(1 + a^) on„2 


(3x3 _ 202.2 _j_ _ 2j log(l -I- x) 


Lg{x) = (2 - 3x^) log(x) 

Lg{x) = (^ + 3^) (^ —11 + a;) _ 2xlog(x) (18) 
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and the rational parts Rx are 


Ra{x) = —— 20a; + 12) 

X 

Rb{x) = (x^ - 10a; + 1) 

X‘^ 

2 9 

Rc{x) = ^ + 2 - a;^ 

X‘^ 

Rg{x) = - + 2. (19) 

X 

In the UV the functions F, G have the asymptotic be¬ 
haviour 

Fi.x) ~ ^ ^ log(a;); G[x) + \ log(x) (20) 

so that the standard UV behaviour is recovered for 

s,So > 1 

J(s)-1 = 1 -p Halog(s/so) 

y(s)-i = i + ia log(s/so) (21) 

as we could predict from the wave function renormaliza¬ 
tion constants Ea. (ll4l) . 

In the opposite limit x —>■ 0 we find G(x) —>■ const, and 
R(x) ^ (1/a;) so that x(0) is finite and J(s) ^ (s/a), 
yielding a hnite gluon propagator A(0) = 8/5am^ as 
predicted by Ea. lfT^ . 

Up to a renormalization factor, the dressing functions 
are invariant for a change of the bare coupling. That is 
more evident if we consider the rescaled functions (aJ) 
and (ay) since then Eas. lfTKll . dlBII loose any explicit de¬ 
pendence on a. In fact, the physical content of the theory 
is inside the universal functions F, G. The actual value of 
takes the role of a natural scale that fixes the physical 
units and can only be determined by a comparison with 
some phisical quantity, as for lattice simulations. That is 
just a consequence of the lack of an energy scale in the 
Lagrangian. However, the present calculation is very pre¬ 
dictive and, up to irrelevant constants, the inverse dress¬ 
ing functions are predicted to take the universal shape of 
the functions F and G: we can write Eqs. ((I5l),((l6l) as 

[ap^A(p)] ^ = F{p^/m^) -\- const. 

[ax(p)]~^ = G{j?In?) F const. (22) 


where the constants depend on normalization, bare cou¬ 
pling and subtraction points. Thus we expect that, up 
to a scaling factor, all lattice data can be put on top of 
the plots of F and G by an additive constant. In Fig.3 
and Fig.4 the lattice data of Ref. [3^ are shown together 
with the plots of the functions F and G. The function F 
has a pronounced minimum at a; ~ 1.62 that is p ~ 0.93 
GeV in the units of the lattice data, thus fixing the scale 
m = 0.73 GeV that is used in the figures. We also show 



Figure 3: The function F(x) c (line) is plotted together 
with the lattice data for the inverse gluon dressing function 
1/J{p) = p^A{p) (points) extracted from the figure of Ref. 

{N = 3, g = 1.02, L=96) and scaled by a renormalization 
factor. The energy scale is set by taking m — 0.73 GeV. 



p (GeV) 


Figure 4: The function G{x) + c (line) is plotted together 
with the lattice data for the inverse ghost dressing function 
l/x{p) (points) extracted from the figure of Ref. [32| {N = 3, 
g — 1.02, L=96) and scaled by a renormalization factor. The 
energy scale is set by taking m = 0.73 GeV. 


the gluon propagator and the ghost dressing function in 
Fig.5 and Fig.6. 

The agreement is very good for a one-loop calculation 
but deviations can be expected when s is far from the 
subtraction point. Thus a slight dependence on the sub¬ 
traction point would be a natural consequence of the one- 
loop approximation. 

Nevertheless, let us pretend that we can ignore such 
limitations and look at observable quantities like the 
running coupling that can be related to phenomenology. 
Assuming that in Landau gauge the ghost-gluon vertex 
is regular}^ and the vertex renormalization constant 
can be set to one in a momentum-subtraction scheme. 
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Figure 5: The gluon propagator A(p) obtained by Eg. (1151) 
(line) is plotted together with the lattice data (points) ex¬ 
tracted from the figure of Ref.j^ {N = 3, g = 1.02, L=96) 
and scaled by a renormalization factor. The energy scale is 
set by taking m = 0.73 GeV. 



(GeV^) 


Figure 7: The running coupling Os = 47ra/(3A) for N = 3 
by Ea. ll23ll . so and Qs(so) are fixed at the maximum (so = 
xm = 1.044, as (so) = 1.9) and m = 0.6 GeV (solid line). The 
points are the lattice data of Ref. IsSl . The broken line is the 
standard one-loop behaviour for A = 0.81 GeV. 



p (GeV) 

Figure 6: The ghost dressing function x(p) (line) is plotted to¬ 
gether with the lattice data (points) extracted from the figure 
of Ref.j^ (N = 3, g = 1.02, L=96) and scaled by a renor¬ 
malization factor. The energy scale is set by taking m = 0.73 
GeV. 

a runnig coupling is usually defined by the renormaliza¬ 
tion group invariant product a(s) = ap J(s)x(s)^ where 
J{so) = x(so) = 1 and op = a(so). Using Eas. dTKl) and 
(USD, the one-loop running coupling can be written as 

" l + a(so)[^(s)-5(so)] ^ ^ 


In the infrared, as shown in Fig.7, the running coupling 
is finite, it does not encounter a Landau pole and vanishes 
in the limit s —1 0 as the power a(s) ~ s. A maximum is 
found at the point where dS'(a;)/da: = 0, which occurs at 
Xm = 1.044. This point does not depend on any parame¬ 
ter, and can be used as an alternative method for fixing a 
physical energy scale if that energy is measured somehow. 
For N = 3, many lattice simulations predict a maximum 
at p Ri 0.6 — 0.7 GeV giving a scale m r; 0.6 GeV that is 
not too far from the value m = 0.73 that was extracted 
from the propagators of Ref.[^. Probably, this small 
difference would be narrowed by the inclusion of higher 
loops in the calculation. 

While in the UV the asymptotic behaviour of the cou¬ 
pling seems to give no information on the scale m, in 
the infrared a phenomenological knowledge of the cou¬ 
pling would fix all the free parameters. Suppose that we 
pinpoint the value of as at its maximum, ag ~ 1.9 ac¬ 
cording to Ref.[^, and let us explore as an exercise the 
behaviour of the coupling when we go back towards the 
UV by Ea. ((2^ . with all the limitations of the one-loop 
approximation. The maximum at p r: 0.6 gives a scale 
m Ri 0.6 GeV. We can set sq at this point Sq = Xm and 
then take as(so) as the maximum coupling. At the same 
point S{xm) = Sm = 3.09. Then, pushing s towards the 
UV and making use of the asymptotic behaviour Eq. (j20ll , 
we can insert S{s) Ri 29/18-1-(11/9) log(s) in Ea. (l23l) and 
write it as a standard one-loop coupling 


where S{x) = F{x) + 2G{x). 

By Ea. (l20l) . for s, sq ^ 1 we find the standard 
UV behaviour a[S'(s) — S'(so)] r: (11q;/9) log(s/so) = 
(11 Vas/127r) log(s/so). In this limit the result does not 
depend on the scale m and we recover the well known 
one-loop running coupling. 




127r 

llVlog (p^/A^) 


where the scale A is defined as 


A = bm exp 


127r 

22Nas{xM) 


(24) 


(25) 
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and the coefficient b satisfies = exp(9S'M/ll—29/22) = 
3.35. Eg. ((2^ provides a direct link between the Landau 
pole of the standard one-loop running coupling at p = A 
and the infrared parameters, namely the mass scale m 
and the maximum value of the coup ling as{xM)- In¬ 
serting the value as(so) ~ 1.9 of Ref. (34| in Ea. d^ we 
obtain A = 0.81 GeV for = 3, which is not too far 
from the value A = 0.7 GeV that is used in the same 
paper for a fit of the lattice data in the UV. Moreover, 
we can extract from Ea. ((25l) the pure theoretical bound 
m > A/&Ri 0.55A. 

Of course, the calculation is just a one-loop approxi¬ 
mation. It generalizes the standard one-loop results to 
the infrared, but mantains the limitations of neglecting 
higher order terms. Moreover, when connecting very 
different scales, the use of renormalization group tech¬ 
niques becomes mandatory for a quantitative description 
and has been shown to be very effective in other massive 
models [si-Iaii. 

In summary, the double expansion has been shown to 
be viable for energies ranging from the UV to the in¬ 
frared, without changing the original Lagrangian, reach¬ 
ing a good agreement with the lattice data from first 
principles. While the full potentialities of the method 
are totally unexplored yet, it might extend the standard 
perturbative appoach to lower energies deep inside the 
non-perturbative sector of QCD. 
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